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1.  STATEMENT  OF  WORK 

-ft" CLOGS' 

The  goal  of  this  research  program  is  an  improved  understanding 
of  the  roles  of  microstructure  and  point  defect  structure  in  deter¬ 
mining  the  mechanical  and  optical  properties  of  infrared  materials 
upon  which  the  development  of  new  materials  and  novel  synthesis 
techniques  can  be  based.  The  work  statement  of  the  program  includes 
the  following  items: 

a.  Determine  the  optical  and  mechanical  properties  of  pure  and 
doped  alkali  halide  specimens  given  various  isothermal  heat  treat¬ 
ments  and  correlate  findings  with  the  microstructures  and  point 
defect  structures; 

b.  Conduct  similar  studies  utilizing  transmission  electron 
microscopy  on  specimens  of  pure  and  doped  alkaline  earth  fluorides; 

c.  Investigate  the  optical  and  mechanical  properties  of  as- 
irradiated  and  annealed  ion- implanted  layers  on  zinc  selenide  and 
correlate  findings  with  microstructure  as  determined  by  trans¬ 
mission  electron  microscopy; 

d.  Study  the  feasibility  of  preparing  a  variety  of  pure  and 
doped  infrared  materials  by  evaporation  driven  liquid  sintering; 

e.  Examine  the  role  of  point  defect  structure  and  microstructure 
in  determining  the  optical  and  mechanical  properties  of  pure  and 

doped  magnesium  oxide  specimens.  Approved  for  public  release; 

distribution  unlimited. 


2.  ACCOMPLISHMENTS 


2.1  Optical  Absorption  Facility  (M.  Bass) 

A  major  facility  for  the  study  of  optical  absorption  has  been 
established.  The  facility  is  shown  schematically  in  Fig.  1.  In 
it,  conventional  spectrophotometry  as  indicated  by  the  Cary  14 
Spectrophotometer  is  supplemented  by  laser  calorimetry.  The  laser 
calorimetry  system  includes  several  CW  lasers,  an  HP  9825  Controller/ 
Calculator  based  data  acquisition,  processing  and  display  system 
and  an  evacuable  calorimeter  chamber.  This  system  can  measure 
with  a  1  W  laser  beam  an  absorption  as  small  as  one  part  in  10^. 

The  lasers  include  a  C02  laser  at  10.6  ym,  an  HF/DF  laser 
at  ^  2.4  -  3.8  ym,  two  Nd:YAG  lasers  (one  at  1.06  and  one  at  1.318  ym) 
a  Kr  ion  laser  and  a  dye 'laser  covering  the  range  from  0.35  to 
1  ym.  The  HF/DF  laser  was  received  late  in  the  first  year  from 
RADC  and  required  an  extensive  overhaul  before  it  could  be 
reliably  operated.  It  was  nearly  ready  for  use  at  the  time  of 
termination  of  this  project. 

2.2  Two  Photon  Absorption  Studies  (M.  Bass,  E.  W.  VanStryland  and 

A.  F.  Stewart) 


Early  in  this  program,  while  preparing  the  laser  calorimetry 
facility,  we  conducted  experiments  concerning  two-photon  absorption. 
This  work  is  described  in  an  "Applied  Physics  Letters"  article. 


which  is  attached  to  this  report  (see  Appendix  1). 
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Fiaure  1  Laser  Absorption  Studies  Laboratory:  Optical  Configuration 


2.3  Materials  (Susan  Allen) 

The  materials  obtained  for  optical  and  mechanical  studies 
are  listed  in  Table  1. 


Composition 

Vendor 

Comments 

CaF2  +1.0% 

Gd 

Raytheon 

Fusion  cast,  coarse-grains 

CaF2  +  0.1% 

Gd3+ 

Optovac 

Crystal 

CaF2  +  1.0% 

Gd3+ 

II  II 

II  II 

CaF2  +1.0% 

Ce3+ 

It  II 

II  II 

SrF2  +  1.0% 

Gd3+ 

It  It 

II  II 

BaF2  +  1.0% 

Gd3+ 

H  n 

II  II 

CaF2  +  Th3+ 

Harshaw 

Crystal  window 

(<  1%  Th3+) 

CaF2  +  Dy 3+ 
KCl  +  Eu2+ 

II  II 

II  II 

II  If 

It  II 

(<  1%  Dy3+) 

KC1  +  Sr2+ 
2+ 

KCl  +  Sr 

CaF2 

NaF 

Raytheon 

II  II 

II  II 

II  II 

•  Crystal  73-82 

"  "  73-63 

Cast  windows 

(5000  PPM  Sr 

Table  1:  Alkali  halides  and  alkaline  earth  halides 
for  optical  and  mechanical  studies 


,2.4  Absorption  Results  by  Laser  Calorimetry  (M.  Bass  and  S.  Allen) 

After  the  nearly  one  and  one  third  years  required  to 
establish  the  optical  absorption  facility,  experiments  were  initiated 
to  measure  the  10.6  and  1.06  pm  absorption  of  samples  of  doped 
and  undoped  alkali  halides,  alkaline  earth  fluorides.  The 
results  of  these  measurements  are  listed  in  Table  2 . 

In  addition  the  2.7  pm  absorption  of  CaF2  (Sample  #  1  in  Table  2  ) 
was  measured  and  found  to  be  4  x  10'^  cm'^. 


Table  2  :  Absorption  of  Alkali  Halides  and  Alkaline  Earth  Halides 


Samples 

C02  10.6 

U 

YAG  1.06  u 

1) 

CaF2  Raytheon 

3.5 

6.5  x  10'5  cm*1 

2) 

NaF 

1.22 

8.5  x  10'5  cm'1 

3) 

1  %  Gd/CaF2 

High* 

6.8  x  10'3  cm'1 

4) 

5000  PPM  Sr++/KCl 
Raytheon  73-82 

4.0  x 

to 

1 

o 

H 

cm  * 

3.6  x  10'4  cm'1 

5) 

Eu++/Kcl  Harshaw  - 

1.13  : 

t  10' 

3  cm'1 

5.48  x  10'4  cm' 

6) 

Th++/CaF2  Harshaw 

High* 

6.9  x  10  4  cm1 

7) 

Dy+++/CaF2  Harshaw 

High* 

.45  cm  1 

8) 

Sr++/KCl 

Raytheon  73-63 

2.45  cm  long 

5.2  x 

10'3 

cm' 1 

3.3  x  10  4  cm"1 

1.65  cm  long 

3.3  x 

10'3 

cm’ 1 

1.5  x  10'4  cm’1 

.96  cm  long 

2.5  x 

10'3 

cm* 1 

.9  x  10  4  cm  1 

* 

Too  large  to  be  measured 


2.5  Vacancy  Tracking  by  Color  Centers  Generated  by  Two  Photon 
Absorption  (H.  Park  and  W.  H.  Steier) 

The  goal  of  this  work  was  to  study  and  understand  the  interaction 
of  the  optical  and  mechanical  properties  of  pure  and  doped  alkali  halide 
crystals.  These  properties  are  correlated  with  point  defects  and  vacancies 
of  these  materials.  As  a  way,  two  photon  excitation  of  electron-hole  pairs 
was  proposed  as  a  mean  of  elucidating  the  point  defect  structure  of  these 
materials.  If  electrons  or  electron-hole  pairs  (i.e.  excitons)  are  excited 
somehow  (e.g.  X-ray,  electron-beam,  or  UV) ,  electrons  will  be  trapped  at 
the  negative  ion  vacancy  site  (holes  at  positive  ion  vacancy  site)  and  pro¬ 
duce  color  centers  (F  or U  centers)  in  the  crystal.  The  density  and  dis¬ 
tribution  of  these  color  centers  should  be  related  to  the  density  and  dis¬ 
tribution  of  the  vacancies. 

During  this  period  we  only  dealt  with  F  centers  because  they  are  re¬ 
latively  easy  to  generate  and  measure.  The  negative  ion  vacancy  density 
should  be  proportional  to  the  F-center  concentration  in  the  crystal  as  a 
first  approximation  and  the  F-Center  concentration  can  be  measured  from 
the  transmittivity  of  crystal  at  the  absorption  band  (F-band)  ^  ^  In  these 
experiments,  a  XeCl  laser  has  been  used  to  generate  electron-hole  pairs  in 
alkali  halide  crystals  by  two  photon  absorption.  Two  photon  energy  of  the 

9 

laser  (3080A)  is  8.06  eV.  which  is  well  within  the  absorption  edge  of  some 

of  alkali  halide  materials.  The  first  absorption  peak  of  NaCl,  KC1  and 

(2) 

KBr  are  7.96  ev,  7.76  and  6.77  eV  respectively  .  These  are  exciton  bands. 
The  band-to-band  transitions  occur  approximately  1  eV  beyond  the  first 
absorption  peak  for  most  ionic  crystals^.  For  all  three  above  crystals. 


color  centers  can  be  easily  obtained  with  XeCl  laser  pulse  energy  less 


than  lOmJ  focused  on  2  x  0.5mm  (pulse  width  -'20r\sec). 

For  these  measurements  KBr  single  crystals  from  Harshaw  Chemical 

Company  were  used.  In  KBr,  bluish  F  centers  were  observed  even  with  a 

5mJ  pulse,  but  on  the  other  hand  these  crystals  are  so  fragile  that  3  or 

4  shots  of  30mJ  pulses  could  easily  damage  the  crystal.  The  absorption 

band  of  the  F  center  in  KBr  was  measured  with  a  spectrophotometer 

°(4) 

and  shown  to  be  centered  around  6300A.  This  wavelength  conveniently  allowed 
the  use  of  the  He-Ne  laser  to  moniter  F-center  concentration  in  KBr. 
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At  t  =  0+  after  0V  laser  pulse,  a  is  13.68  cm-^  from  the  data  which  gives 
an  F  center  concentration  of  5.6  x  10  cm  J  The  decay 

time  of  color  center  is  -'10  sec  from  this  experiment. 

In  the  next  experiment,  we  attempted  to  determine  if  the  F-center  con¬ 
centration  can  be  related  to  the  number  of  vacancy  sites  in  the  crystal. 

In  these  preliminary  experiments  we  attempted  to  vary  the  vacancy 
concentration  by  heating  and  quenching  the  samples.  It  is  known  that  one 
can  generate  additional  vacancy  sites  in  a  crystal  by  heating  it  close  to 
the  melting  point,  and  that  these  vacancies  can  be  frozen  in  by  rapidly 
quenching  the  crystal.  In  this  way  one  can  increase  the  vacancy  concen¬ 
tration  in  a  crystal however  the  excessive  vacancies  will,  with  time, 
diffuse  to  the  surface  of  the  crystal ^  .  In  this  experiment  we  measure 
the  transmittivity  of  the  crystal  before  and  after  quenching  at  t  -  0+ 
after  exposure  to  the  same  UV  energy 

The  KBr  crystal  was  held  at  650°C  for  20  minutes  in  a  oven  and  then 
quenched  to  room  temperature  in  oil.  The  transmission  was  measured  at 
t  =  0+  as  a  function  of  time  after  the  quench.  The  results  are  shown  in 
Figure  4.  Each  circle  is  the  average  value  of  5  measurement  at  different 
spots  of  the  crystal. 

From  the  data  we  can  conclude  that  up  to  1/2  hour  after  quenching  the 
transmittivity  (or  F  center  concentration)  is  different  from  that  before  the 
heating  and  quenching  but  after  one  hour  the  transmittivity  has  returned 
to  its  original  value.  The  transmittivity  change  is  ~3%  which  is  marginally 
greater  than  our  experimental  error.  In  another  experiment  with  another 
KBr  crystal,  the  transmittivity  decreased  by  *'■5%  up  to  1-1/2  hours  after 


the  quench. 


after  quenching) 


TIME  AFTER  QUENCHING  (in  hours) 


f 


We  can  tentatively  interpret  our  results  to  confirm  the  relationship 
between  vacancy  density  and  color  center  density,  however  other  factors 
must  also  be  considered.  For  example,  it  is  known  that  in  crucible  grown 
KI  the  color  center  density  is  a  function  of  previous  irradiation  and 
plastic  deformation,  however,  other  work  has  shown  this  not  to  be  important 
in  KBr.(s) 

These  preliminary  experiments  appear  to  be  successful  but  considerable 
more  work  remains  to  show  that  the  color  center  density  is  a  valid  measure 


of  vacancy  density. 
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2.6  Optical  and  Mechanical  Properties  of  Alkali  Halides 

and  Alkaline  Earth  Fluorides  (S.  M.  Copley  and  J.  Savage) 

Because  so  little  research  had  been  done  on  the  effect  of 
doping  on  the  mechanical  properties  of  the  alkaline  earth  fluorides, 
it  was  decided  to  initiate  work  on  this  type  of  material.  Alkali 
halide  specimens  were  also  obtained,  see  Table  1,  however,  only 
optical  measurements  had  been  carried  out  on  these  specimens 
at  the  time  this  program  terminated. 

The  goal  of  this  research  was  to  determine  by  hardness 
and  absorption  measurements  (both  I.R.  and  U.V.)  whether  thermal 
treatment  could  be  used  to  harden  alkaline  earth  fluorides  through 
the  formation  of  point  defect  complexes  without  producing  an 
increase  in  I.R.  absorption.  The  U.V.  absorption  measurements 
were  to  be  carried  out  -as  a  means  of  detecting  defect  complex 
formation. 

Experiments  were  initiated  on  CaF2  +  1.0%Gd^+  deformed 
bend  specimens  provided  by  Raytheon.  The  specimens  had  been 
originally  been  cut  from  fusion  cast  ingots  and  were  very 
coarse-grain  polycrystals. 

An  initial  experiment  was  carried  out  on  an  as-received 
specimen  to  determine  the  amount  of  scatter  that  could  be 
expected  in  a  series  of  hardness  measurements.  A  traverse 
was  made  along  one  face  of  a  specimen.  Indentations  were  made 
with  the  long  dimension  either  parallel  to  or  perpendicular  to 
the  long  dimension  of  the  specimen.  The  measurements  were  made 
over  a  period  of  40  days  to  detect  any  change  with  time.  The 
best  estimate  of  the  population  mean  and  standard  deviation 
was  214.4  and  a  =  5.6,  respectively,  for  the  transverse  measure- 


merits  and  \i  =  201.6  and  o  *  1.4  respectively  for  the  longitudinal 
measurements.  It  was  concluded  from  these  measurements  that 
the  scatter  was  small  enough  so  that  hardness  could  be  used  to 
measure  any  significant  change  in  yield  stress  and  that  neither 
location  along  the  specimen  nor  time  produced  any  systematic 
variation  in  hardness . 

To  investigate  the  effect  of  thermal  treatment,  specimens 
were  annealed  at  625 °C  to  dissociate  all  point  defect  complexes, 
quenched  and  then  aged.  It  was  hoped  that  the  quench  would  be 
rapid  enough  to  preserve  the  point  defects  in  a  dissociated 
state  but  that  the  defects  would  form  complexes  upon  aging 
resulting  in  an  increase  in  hardness.  It  was  found,  however, 
that  direct  quenching  into  oil  was  not  possible  because  of 
extensive  cracking  due  to  thermal  stresses.  Consequently, 
experiments  were  carried  out  in  which  the  specimen  was  sealed 
in  a  fused  quartz  capsule  during  oil  quenching. 

Table  3  shows  the  results  of  several  quenching  experiments. 
Each  hardness  value  is  the  average  of  five  independent  determina¬ 
tions.  Although  the  results  are  not  entirely  consistent,  it 
appears  that  annealing  and  quenching  increases  hardness.  This 
suggests  that  the  quench  was  not  sufficiently  rapid  to  suppress 
the  formation  of  hardening  point  defect  complexes. 

Although  it  was  originally  intended  to  age  the  quenched 
specimens  at  a  temperature  greater  than  room  temperature,  it 
was  discovered  that  changes  in  hardness  occurred  during  room 
temperature  aging.  Typically,  the  hardness  decreased  initially 


Sample 

As-Received  (KHN) 

Quenched 

(KHN) 

L 

T 

L 

T 

9  a 

242 

215 

262 

(+20) 

265 

(+50) 

13  a 

226 

242 

222 

(-4) 

195 

(-47) 

14  a 

227 

232 

216 

(-11) 

240 

(+8) 

15  a 

204 

204 

229 

(+25) 

238 

(+34) 

16  a 

216 

214 

223 

(+7) 

218 

(+4) 

17  a 

208 

227 

225 

(+17) 

226 

(-1) 

Table  3; 

Knoop  hardness 

of  various 

specimens 

in 

as-received  and  quenched 

conditions. 

The  letter  L 

indicates  that  long  axis  of  indent  was  parallel  to 
the  long  axis  of  the  specimen;  T  indicates  that  the 
long  axis  of  the  indent  was  perpendicular  to  the  long 
axis  of  the  specimen. 


and  then  increased  passing  through  a  maximum  at  about  50  hours. 

The  scatter  in  results  of  the  quenching  experiments  as  well 
as  the  occurrence  of  "hardness  reversals"  in  the  aging  curves 
suggests  that  the  Raytheon  samples  may  include  several  grains. 

In  this  case,  deviations  from  normal  behavior  could  be  attributed 
to  an  orientation  effect  or  to  the  nearness  of  the  indentation 
to  a  grain  boundary. 

Preliminary  examination  of  the  uv  spectra  of  the 

Raytheon  specimens  after  quenching  and  agin  produced  no  evidence 

for  defect  complex  formation.  Thus  while  we  have  shown  that 

3+ 

thermal  treatment  may  increase  the  hardness  of  Gd  doped  CaF2 

at  room  temperature,  at  least  temporarily,  we  have  not  been 
able  to  link  the  hardness  increase  with  the  formation  of  any 
specific  defect. 


2.7  Optical  and  Mechanical  Properties  of  ZnSe  (W.  G.  Spitzer  and  J.  S.  Ko 


A  manuscript  describing  the  results  of  this  work  has  been 
prepared  for  submission  to  the  Journal  of  Applied  Physics  •  A 
copy  of  this  manuscript  is  attached  to  this  report  in  Appendix  B. 

2.8  Evaporation  Driven  Liquid  Sintering 

In  this  project,  the  feasibility  of  applying  a  new  process, 
evaporation  driven  liquid  sintering  (EDLS)  to  the  sintering  of  MgF2 
was  investigated.  Sodium  chloride  was  selected  as  the  volatile 
liquid  sintering  aid.  Magnesium  fluoride  powder  compacts  containing 
NaCl  were  heated  in  a  controlled  atmosphere  to  suppress  evaporation 
of  NaCl.  After  sintering  has  occurred,  the  NaCl  must  be  evaporated 
from  the  compact. 

It  has  been  shown  that  for  the  evaporation  driven  liquid 
sintering  process  to  be  successfully  applied,  the  volatile  liquid 
must  completely  wet  the  the  powder  being  sintered  and  the  powder 
being  sintered  must  have  limited  solubility  in  the  volatile  liquid. 
Neither  of  these  conditions  were  satisfied  by  the  NaCl.  Before 
another  sintering  aid  could  be  found,  the  student  working  on 
this  project  left  the  University. 
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MgO  is  a  mixed  conductor  the  properties  of  which  are  strongly 
affected  by  aliovalent  dopants  such  as  Fe,  Mn  and  Li.  The  con¬ 
tributions  by  ionic  and  electronic  components  to  the  current  can  be 
determined  by  combining  measurements  of  total  conductivity  with 
transference  numbers  t^  and  tgl  *  1  -  t^  determined  from  emf 
measurements  on  cells  of  the  type: 

Pt,02  (p02jI)/Mg0/Pt,02(p02)n 

the  emf  of  which  is  given  by: 
c  RT  11 

E  =  ?F  /  -t.  d  In  Po, 

I  1  2 

Transference  numbers  at  (^2^11  are  therefore  found  from  the  slope  of  a 
curve  versus  (P°2^n  wit^  (P02)  I  constant: 

CVn  -  £  (3E/31n(P02)n) 

1  u  KT  2  11  (p02)I 

and  °i  =  %  ’  ael  =  a  {1'ti) 

A  high  temperature  furnace  with  super  Kanthal  heating  elements 
and  Zircon  insulation was  built,  a  gas  train  for  mixing  gases 
required  for  establishing  oxygen  pressures  between  10  10  and  1  atm 
was  constructed,  and  measuring  equipment  was  assembled.  Owing  to 
personnel  problems,  no  measurements  were  done  before  the  contract 
expired. 
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Laser  calorimetric  measurement  of  two-photon  absorption 

Michael  Bass,  Eric  W.  Van  Stryland,  and  A.  F.  Stewart 
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A  laser  calorimeter  has  been  used  to  measure  two*photon  absorption  in  solids.  This  new 
technique  provides  greater  sensitivity  than  is  possible  in  nonlinear  transmission 
measurements.  The  two-photon  absorption  coefficients  of  CdTe  and  CdSe  at  1.06  pm 
using  ~  16-nsec  (FWHM)  pulses  are  0.13j;0.04  and  0.050^-0.014  cm/MW, 
respectively. 


PACS  numbers:  42.10.Ke,  78.20.Dj,  42.60. -v 

The  laser  calorimeter  is  a  device  that  is  widely  used  to 
measure  small  linear  absorptions  at  wavelengths  where  laser 
sources  are  available.  Such  an  instrument  can  measure  ab¬ 
sorptions  which  are  less  than  one  part  in  105  using  a  1-W 
laser.'  The  calorimeter  measures  the  total  absorption  due  to 
all  processes  including  the  intensity-dependent  phenomenon 
of  two-photon  absorption  (TP A).  Thus,  the  experimental 
technique  described  in  this  paper  is  based  on  a  measurement 
of  the  total  absorption  as  a  function  of  the  laser  intensity. 
Calorimetric  measurements  of  the  total  absorption  in  CdTe 
and  CdSe  as  functions  of  1.06-ytm  laser  intensity  were  used 
to  obtain  both  linear  and  TPA  coefficients  of  these  materials. 
The  results  can  be  understood  by  using  a  simple  model  for 
the  attenuation  with  distance  in  the  presence  of  TPA  and  by 
properly  accounting  for  multiple  reflections  in  the  sample. 

In  1931,  M.  Goppert-Mayer  presented  a  theoretical  de¬ 
scription  of  TPA  in  which  the  transition  probability  was 
shown  to  be  proportional  to  the  square  of  the  light  intensity.' 
Thus,  a  material  having  both  linear  and  two-photon  absorp¬ 
tion  will  transmit  light  according  to  the  expression 

^=fx,yj,t )  =  -[a +01  (x,yj,t )]/  {x,yjt,t ).  ( 1 ) 

dZ, 

Here,  I  (x.ya,t )  is  the  light  intensity  in  MW  /cm2  and  is  a 
function  of  the  space  coordinates  and  time,  z  is  the  direction 
of  propagation,  a  is  the  linear  absorption  coefficient  in  cm'1, 
and  0  is  the  two-photon  absorption  coefficient  in  cm/MW. 
In  many  materials  of  interest  a  will  be  small,  a~10"!  cm'1, 
and/?~10':  cm/MW  according  to  the  transmission  mea¬ 
surements  at  high  intensities. '  Thus,  a  modest  intensity  of  l 
MW/cnr  will  result  in  a  two-photon  contribution  to  the  to¬ 
tal  absorption  which  is  equal  to  the  linear  term  [see  Eq.  (I)]. 
On  the  other  hand,  this  intensity  will  result  in  a  change  in 
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FIG  1  Experimental  apparatus  for  TPA  measurements  using  laser 
calorimetry 


transmission  for  samples  shorter  than  1  cm  of  less  than  1%, 
which  is  too  small  to  be  measured  accurately.  A  comparison 
of  the  two  techniques  can  be  made  by  noting  that  in  a  trans¬ 
mission  experiment  one  measures  a  small  change  (the  reduc¬ 
tion  of  transmission  due  to  TPA)  in  a  large  quantity  (trans¬ 
mitted  intensities),  whereas  using  calorimetry  one  measures 
a  large  change  (the  added  absorption  due  to  TPA)  in  a  small 
quantity  (total  absorption). 

The  experimental  procedure  is  to  measure  the  total  ab¬ 
sorption  as  a  function  of  intensity  and  evaluate  the  data  ac¬ 
cording  to  Eq.  (1),  taking  into  account  reflections  at  the  rear 
surface.  A  flashlamp-pumped  ^-switched,  TEM«,  mode 
Nd  :  YAG  laser  was  the  source  for  the  experiments,  as 
shown  in  Fig.  1 .  The  laser  was  operated  at  —  5  Hz  with  an 
output  of  —  1  mJ/pulse.*  The  pulse  temporal  waveform  was 
a  Gaussian  function  with  full  width  at  i/e  points  in  intensity 
of  19  nsec.  A  beam  scan  was  used  to  determine  the  spatial 
profile  which  was  Gaussian,  having  a  1/e  full  width  in  inten¬ 
sity  of  1.2  mm  at  the  focusing  lens.  An  8-cm  focal  length  lens 
was  used  to  increase  the  intensity  at  the  samples,  which  were 
placed  4  cm  behind  the  focus.  The  attenuators  and  optics  in 
front  of  the  sample  limited  the  maximum  peak-on-axis  inten¬ 
sity  at  the  sample  to  ~  30  MW  /cm2,  with  an  average  power 
of  approximately  5  mW.  The  energy  per  pulse  was  measured 
by  averaging  the  output  of  over  100  pulses  with  a  pyroelec¬ 
tric  joulemeter.  The  accuracy  of  the  TPA  measurements  is 
sensitive  to  pulse-to-pulse  energy  variations,  which 
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FIG  2.  Intensity-dependent  absorption  in  CdSe:  high-mtenstty  data  (1 _ Nt; 

low-intensity  tew  I  data  (  ):  theory  <-) 
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FIG.  3.  Intensity-dependent  absorption  in  CdTe:  high-intensity  data  (Us); 
low-intensity  (ew)  data  ( X );  theory  (-). 


must  therefore  be  kept  as  low  as  possible  (  ~  5%).  The  beam 
intensity  at  the  sample  was  adjusted  using  a  two-polarizer 
attenuator  described  in  Refs.  5  and  6. 

A  thermocouple  attached  to  the  otherwise  thermally 
isolated  sample  recorded  the  sample’s  temperature  change 
due  to  the  absorption  of  the  laser  light.  When  amplified  and 
recorded,  the  thermocouple  voltage  appeared  as  shown  on 
the  chart  recorder  in  Fig.  1.  While  the  laser  is  on  the  tem¬ 
perature  rises.  The  gentle  slopes  while  the  laser  is  off  are  due 
ter  small  thermal  losses  that  are  seen  with  thin  samples.  The 
time  constant  of  the  calorimeter  was  long  enough  so  that  a  5- 
Hz  laser  beam  caused  no  uneveness  in  the  recorded  tempera¬ 
ture  change.  Thus,  the  calorimeter  responded  to  the  total 
average  power  absorbed.  The  calorimeter  was  of  convention¬ 
al  design  and,  when  evacuated,  had  a  measured  sensitivity  to 
absorption  of  one  part  in  105  with  a  1-W  source. 

Data  such  as  in  Fig.  1  was  obtained  at  several  laser 
intensities  for  room-temperature  samples  of  CdTe  and  CdSe. 
The  samples  were  single  crystals,  and  the  CdSe  was  oriented 
with  its  optic  axis  normal  to  the  entrance  (exit)  surface.  This 
orientation  eliminated  the  possibility  of  second-harmonic 
generation  at  0.53/jm  which  could  be  directly  absorbed.  The 
CdTe  sample  had  its  optic  axis  approximately  22s  away  from 
the  normal  to  the  entrance  surface.  Figures  2  and  3  show  the 
intensity  dependence  of  the  total  absorption  of  these  two 
materials  measured  calorimetrically.  For  a  small  total  ab¬ 
sorption  [i.e.,  (a  +01  )z<  1  ],  Eq.  ( l )  predicts  a  straight  line  for 
absorption  versus  intensity.  The  vertical  intercept  of  this  line 
gives  a  and  its  slope  determines  0.  In  these  experiments 
which  employ  a  Gaussian  beam  profile  and  a  Gaussian  tem¬ 
poral  shape,  the  TPA  coefficient  would  be  2%/  2  times  the 
slope  in  such  a  plot.  'Hie  data  fits  a  straight  line  quite  well, 
but  the  total  absorption  at  the  highest  intensities  used  for 
CdTe  was  over  20%.  Thus,  Eq.  ( 1)  must  be  solved  exactly. 
This  was  done  previously,  including  the  effects  of  Gaussian 
spatial  and  temporal  profiles,  but  the  effects  of  multiple  re¬ 
flections  were  neglected.'  For  high-index  materials,  the  in¬ 
tensity  of  the  reflected  beam  from  the  rear  surface  can  sig¬ 
nificantly  affect  the  computed  absorption  coefficients.  Since 
TPA  results  from  the  imaginary  part  of  the  third-order  opti¬ 


cal  nonlinear  susceptibility  ^(J\  the  field  equations  for  both 
the  incident  and  reflected  beams  were  used  to  obtain  the 
following  coupled  differential  equations  for  the  intensities: 

^-=-{a+0(I+IR)]I, 

dz 

—  =  +  [a  +0(1 +/*)]/*, 

where  /  is  the  intensity  of  the  incoming  beam  and  Ig  is  the 
intensity  of  the  reflected  beam.  Subsequent  reflections  have 
been  neglected  because  they  contribute  still  less  to  the  total 
absorption.  Noting  that  the  product  I  (z\IR(z)  is  a  constant 
everywhere  in  the  sample,  these  equations  can  be  solved  in 
closed  form  with  the  constant  as  a  parameter.  A  rapidly  con¬ 
verging  numerical  iterative  procedure  for  determining  the 
constant  from  the  boundary  conditions  was  used  to  deter¬ 
mine  the  total  absorption  integrated  over  the  spatial  and 
temporal  beam  profiles. 

The  theoretical  fits  along  with  the  data  are  shown  in 
Figs.  2  and  3.  These  fits  were  used  to  obtain  the  values  for/? 
listed  in  Table  I.  The  values  given  fora  were  obtained  using  a 
cw  Nd  :  YAG  laser  and  the  usual  calorimetric  technique. 
These  measured  values  of  a  are  quite  accurate  and  as  such 
were  used  in  the  calculations  of/9.  This  means  that  the  values 
of/9  given  in  Table  I  are  the  result  of  a  one-parameter  fit.  The 
deviation  of  the  fit  from  the  data  at  high  intensities,  particu¬ 
larly  for  CdSe,  may  be  due  to  a  contribution  from  direct 
absorption  by  two-photon-created  free  carriers,  as  discussed 
in  Refs.  3  and  7.  However,  at  this  time  the  data  is  insufficient 
and  the  error  bars  too  large  to  make  any  conclusions  con¬ 
cerning  the  possible  role  of  free  carriers.  Work  is  continuing 
in  this  area  seeking  a  predicted  pulse  width  dependence  to 
TPA  if  free-carrier  absorption  is  significant. 

The  experimental  error  bars  shown  in  Figs.  2  and  3  were 
determined  from  scatter  in  the  data  and  uncertainties  in  de¬ 
termining  the  heating  rate  from  curves  such  as  sketched  on 
the  chart  recorder  in  Fig.  1 .  The  laser-intensity  fluctuations, 
the  accuracy  of  the  laser-spot-size  measurement,  and  tempo¬ 
ral  profile  measurement  along  with  the  uncertainties  men¬ 
tioned  above  combine  to  give  the  quoted  accuracies  in  Table 
I.  The  accuracy  of  TPA  measurements  using  laser  calori¬ 
metry  has  the  potential  for  considerable  improvement.  The 
initial  work  presented  here  has,  however,  demonstrated  the 
feasibility  of  the  technique. 

The  inclusion  of  the  reflected  beam  in  the  analysis  of 
TPA  measurements  is  essential  in  high-index  materials, 
where  it  has  been  ignored  previously,'  but  may  be  eliminated 
in  future  experiments  (thus  greatly  simplifying  the  calcula¬ 
tions)  by  using  AR-coated  samples.  It  should  be  pointed  out 


TABLE  I.  Linear  absorption  coefficient  <7  and  two-photon  absorption  coef¬ 
ficient  (3  at  I  06  /ini. 


Material 

nr  (cm  ) 

/}  (cm/MW) 

CdSe 

0  062  —0.0b 

0050-0014 

CdTe 

0  335—0.03 

0  130-0  036 
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that  ih  certain  materials  the  two-photon  excited  level  can 
decay  through  a  combination  of  radiative  and  nonradiative 
processes.  Then,  part  of  the  total  TPA  may  not  register  as 
heat  generated  in  the  sample.  This  means  that  a  simulta¬ 
neous  transmission  measurement  may  make  it  possible  to 
determine  the  quantum  efficiency  of  such  radiative  transi¬ 
tions.  This  type  of  quantum  efficiency  study  and  measure¬ 
ments  of  TPA  using  two  different  lasers  are  in  progress. 
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I  Introduction 


There  has  been  a  long-standing  interest  in  the  optical  properties 
of  ZnSe,  particularly  the  luminescent  characteristics  and  their  relation¬ 
ships  to  impurities  and  native  defects.  Some  defect  models  have  been  sug¬ 
gested,  but  the  microscopic  structure  is  not  yet  completely  clear.  Recent¬ 
ly  the  importance  of  ZnSe  as  an  infrared  window  material  for  high  power 
laser  applications  demonstrated  the  need  more  detailed  studies  of  the 
micros true tural  features  of  the  defects  and  their  relationships  to  the 
infrared  absorption. 

It  is  now  well^  established  that  if  the  mass  of  a  substitutional  dopant 
is  sufficiently  small  compared  to  those  of  the  host  atoms,  such  as  A1 
or  Li  in  ZnSe,  then  infrared  absorption  measurements  of  the  localized 
vibrational  modes  (LVM)  are  a  very  useful  diagnostic  tool  in  helping 
to  understand  the  nature  of  the  defect.  An  analysis  of  the  infrared 
spectra  can  give  important  information  about  the  microstructure  and  sym¬ 
metry  properties  of  the  defects  and  impurity  atom  positions  in  the  crystal 
lattice.  Lithium,  with  its  light  mass  and  fast  diffusion  rate  has  played 

a  particularly  important  role  in  LVM  studies  of  several  defects.  Li  dif- 
2 

fused  into  GaAs  can  be  a  n-type  or  a  p-type  impurity,  and  can  sit  at  an 

interstitial  site,  at  a  substitutional  site,  or  in  the  form  of  complexes, 
physical 

Another^property  which  is  impurity  dependent  is  the  hardness  of  materials. 

3 

It  is  known  that  hardening  of  Si-doped  GaAs  is  produced  by  the  interaction 
of  the  solute-vacancy  pairs  with  the  moving  dislocations.  Because  of  the 
IR  window  application,  there  is  considerable  interest  in  hardening  mechanisms 
in  ZnSe.  If  one  labels  the  defects  by  their  LVM,  then  one  can  determine 
whether  the  hardness  of  ZnSe  is  a  function  of  impurity  concentration  and 


if  so,  how  Che  hardness  depends  upon  Che  concencracion  of  specific 
impuriCy-re laCed  defecCs. 

This  is  a  reporc  of  Che  progress  in  a  sCudy  of  Che  opcical  and 
mechanical  properdes  of  ZnSe,  and  ic  includes: 

1)  Localized  vibracional  mode  absorpdon  and  hardness  sCudies  of  Al- 
doped  ZnSe.  BoCh  properties  are  measured  as  functions  of  Impurity 
concentration  and  annealing  conditions. 

2)  A  LVM  study  of  lithium  diffused  ZnSe  which  is  initially  undoped 
or  doped  with  n-type  impurities  (Al,  Ga). 

II  Experimental  Method 

The  ZnSe  ingots  used  in  this  study  were  purchased  from  Eagle- 
Picher  Co.  The  crystals  were  grown  under  an  inert  gas  pressure  by  using 
a  modified  Bridgeman  technique.  The  ingots  are  polycrystalline  with  large 
grain  size  with  the  size  depending  somewhat  upon  the  dopant  and  its  con¬ 
centration.  The  dopants  and  their  concentrations  in  the  material  we  have 

4 

used  are  as  follows  (values  given  by  the  supplier  ) : 

(1)  undoped  ZnSe 

(2)  Al  doped  ZnSe:  100  ppm  and  300  ppm 

(3)  Li  doped  ZnSe:  100  ppm 

(4)  Ga  doped  ZnSe:  100  ppm  and  3000  ppm 

(5)  In  doped  ZnSe:  40  ppm  and  400  ppm. 

The  group  III  elements  Al,  Ga  and  In  are  expected  to  be  n-type 
dopants  In  ZnSe.  The  impurities  are  known  to  be  substitutional  on  the 
Zn  sublattice  and  form  shallow-state,  effective  mass-type  donor^.  Very 
little  is  known  about  the  nature  of  Li  as  an  impurity  in  ZnSe.  An  ac¬ 
ceptor  state  associated  with  Li^  has  been  identified^,  however  the  be¬ 
havior  of  Li  in  ZnSe  is  likely  to  be  complicated  and  involve  more  than 


one  type  of  defect  species.  The  as-grovn  crystals  are  highly  self- 
compensated  and  have  no  detectable  free  carrier  absorption. 

We  used  the  most  heavily  41 -doped  ZnSe  (300  ppm)  in  a  number  of 

annealing  experiments  in  which  the  Se  pressure  is  controlled.  The  samples 
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were  cut  to  a  size  of  about  12  x  5  x  0.3  mm  ,  lapped,  and  then  mechanical¬ 
ly  polished.  After  subsequent  cleaning,  the  samples  were  sealed  in  evacu- 
lated  quartz  ampoules  containing  Se  shot.  The  ampoule  was  placed  in  a 
two  zone  furnace  with  the  crystal  at  one  end  of  the  ampoule  which  was  in 
the  center  of  the  main  furnace.  The  Se  pool  was  at  the  other  the  of  the 
ampoule  which  was  placed  in  the  center  of  the  2nd  furnace.  The  temperature 
of  the  main  furnace  was  set  at  900  K.  The  Se  partial  pressure  was  con¬ 
trolled  by  adjusting  the  temperature  of  the  secondary  furnace.  After  an¬ 
nealing,  the  sample  was  lapped  again  on  both  sides,  reduced  to  a  suitable 
thickness,  and  chemically  polished. 

Infrared  absorption  measurements  were  made  at  liquid  nitrogen  temp¬ 
erature  by  using  a  Perkin-Elmer  model  210  single-beam  grating  spectro¬ 
meter.  The  frequency  range  was  from  v  'V  500  cm  ^  to  'v  300  cm  ^  and  the 
typical  thickness  of  the  samples  measured  was  8  mils  (0.2  mm).  The  absorp¬ 
tion  coefficient  a  was  calculated  from  the  transmission,  T,  and  the  reflect¬ 
ivity,  R,  by  using  the  usual  expression  for  a  solid  sample. 

The  lithium  diffusion  was  carried  out  from  an  alloyed  layer  by  the 
same  technique  as  used  previously  for  GaAs.  In  brief,  thick  samples  were 
surface  alloyed  with  a  mineral-oil  suspension  of  fine  Li  powder  at  500°C 
in  an  argon  atmosphere.  A  piece  of  undoped,  highly  polycrystalline  ZnSe 
which  was  grown  by  Raytheon  Corp.  was  used  to  provide  the  alloy  pieces 
used  in  the  diffusion  process.  The  diffusion  was  carried  out  by  sandwich¬ 
ing  the  samples  between  the  alloyed  pieces.  The  diffusion  temperature 


Q  O 

was  between  600  C  and  800  C.  For  most  samples,  the  time  and  temperature 
used  for  the  diffusion  is  essentially  the  same,  independent  of  dopant  and 
concentration,  about  three  to  four  hours  at  700° C.  As  will  be  discussed 
there  is  evidence  that  one  Li  defect  species  does  increase  with  longer 
diffusion  times  and  therefore  a  Li  saturation  condition  was  not  achieved. 

O 

Samples  diffused  at  800  C  for  longer  than  two  hours  are  heavily  damaged. 
Often  the  central  piece  in  the  sandwich  was  fused  to  the  alloy  piece  or 
the  pieces  became  very  brittle.  After  diffusion  annealing,  the  samples 
were  lapped  and  polished  in  the  usual  manner  and  the  infrared  transmission 
was  measured. 

Hardness  measurements  were  carried  out  for  some  samples  with  a  "Tukon 
Tester"  by  using  a  Knoop  diamond  indentor  loaded  with  a  50  gm  weight. 
Samples  with  chemically  polished  surfaces  were  used.  Indentations  were 
made  on  the  central  area  of  each  grain  to  minimize  the  grain  boundary 
effect.  For  samples  which  were  annealed,  the  hardness  was  measured 
before  and  after  repolishing  the  surface  and  no  change  was  observed. 


III.  Results 


A)  Annealing  of  ZnSe;Al 

The  Infrared  absorption  spectra  of  ZnSe:Al  was  first  studied  by 
Mitsuishi  et  al£  Six  local  inodes  were  observed  due  to  the  presence  of 
A1  as  shown  in  Fig.  1  for  a  sample  having  an  A1  concentration  of  ^  300ppm. 

The  LVM  studies  by  Mitsuishi  et  al.  of  absorption  as  a  function  of  A1 
concentration  and  of  ZnSetAl  annealing  in  a  Se-rich  atmosphere  indicated 
a  different  behavior  for  the  intensity  of  one  of  the  bands  which  has  its 
absorption  peak  at  v  ■  359  cm  \  This  band  was  attributed  to  the  AlZn 
defect,  where  means  X  occupying  a  y  atom  site.  The  source  of  the 
other  five  Al-related  bands  remains  unresolved.  Later  experiments  study¬ 
ing  ZnSe:Al  counterdoped  with  Cu,  Ag,  and  Au  also  confirmed  the  different 
character  of  the  359  cm  ^  band  and  supported  the  AlZn  idendification. 

Mitsuishi  et  al.  attempted  to  obtain  more  information  about  the  Al- 
related  bands  by  studying  the  annealing  ZnSe:Al  in  a  Se-rich  atmosphere. 
Mitsuishi  et  al.  annealed  ZnSe:Al  at  900°  C  for  120  hours  with  a  Se  pressure 
of  2  atms.  All  of  the  Al-related  bands  disappeared  except  the  359  cm 
band.  They  attributed  this  result  to  the  outdif fusion  of  the  Al. 

In  the  present  study  the  annealing  experiments  were  repeated  under 
the  same  conditions  except  shorter  annealing  times  were  used.  Figure  1 
shows  the  absorption  spectrum  of  ZnSerAl  (300ppm)  before  annealing,  after 
aging  for  40  hours  and  80  hours  at  900°C  with  a  Se  pressure  of  2  atms. 

The  359  cm  *  band  increases  initially  during  the  first  40  hours  of  an¬ 
nealing  and  then  remains  unchanged  through  the  remaining  heat  treatment. 

The  five  unknown  bands  all  decrease  with  annealing  time.  Additional 
annealing  experiments  under  different  annealing  conditions  show  the  same 
pattern  with  the  amount  of  change  depending  on  the  temperature  and  Se  pressure. 
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B)  Undoped,  Li  diffused  ZnSe 

Infrared  studies  of  ZnSe  which  were  doped  with  Li  during  crystal 
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growth  have  been  reported  by  Mitsuishi  et  al.  Figure  2 -(a)  and  2 -(b) 
are  approximate  line  spectra  indicating  the  frequencies  and  strengths 
of  the  absorption  bands  observed  by  them.  They  identified 
four  ^Li  bands  and  two  ^Li  bands  as  indicated  by  the  arrows.  They  pro¬ 
posed  that  the  other  bands  between  V  *  450  cm  ^  and  V  *  345  cm  ^  are 
the  result  of  defect-induced  enhancement  of  the  usual  two  phonon  lattice 
bands  and  therefore  did  not  consider  them  further.  They  assigned  the 
bands  at  v  •  336  and  328  cm  ^  to  ^Li  defect  modes  and  the  corresponding 
bands  at  V  *  350  and  341  cm  ^  to  ^Li  modes.  The  isotope  frequency  ratio 
is  1.04  which  is  small  compared  with  1.072  for  Li  in  GaAs.  This  observa¬ 
tion  would  indicate  that  the  local  modes  of  Li  in  ZnSe  consist  of  vibra¬ 
tions  of  polyatomic  complexes  containing  Li  and  a  significant  portion 
of  the  kinetic  energy  of  each  mode  is  not  in  Li  motion.  This  result  is 
unlike  that  for  the  modes  in  Li-diffused  GaAs  which  consist  primarily 
of  Li  vibrations. 

We  used  the  technique  described  in  the  experimental  section  to  dif¬ 
fuse  Li  into  ZnSe.  In  the  Eagle-Dicher  Co.  material, which  was  diffused 
at  700°C  for  4  hours,  we  observed  six  new  ^Li  bands  with  peak  absorption 
frequencies  of  V  =  420,  412,  389,  359,  350  and  343  cm  ^  as  well  as  their 
corresponding  bands  for  ^Li  at  V  =  451,  441,  419,  385,  375  and  367  cm 
(see  Figures  2  (c)  and  2(d)).  The  isotope  frequency  ratio  is  1.072  which  is 
the  same  as  that  of  Li-diffused  GaAs  and  close  to  1.080  for  a  purely  Li  har¬ 
monic  oscillator.  This  shift  indicates  that  the  local  modes  which  we  have 
observed  consist  primarily  of  Li  vibrations.  The  \i  band  at  v  *  420  cm 
as  well  as  the  6Li  bands  at  v  »  441  and  419  cm-1  fall  on  the  strong  two-phonon 
lattice  absorption  bands  but  the  dramatic  increase  in  the  absorption  co- 


efficients  and  the  presence  of  the  isotopic  counterparts  make  it  clear 
that  these  are  LVM  bands  and  not  an  enhancement  of  the  ZnSe  two-phonon  abr 
sorption  process.  We  also  observed  two  additional  weak  Li  bands  at 
v  ■  339  and  326  cm  The  corresponding  \i  bands  are  not  observed, 
probably  because  the  frequencies  are  lower  than  320  cm  ^  where  the  lat¬ 
tice  absorption  increases  rapidly  with  decreasing  frequency  making  ob¬ 
servation  of  LVM  absorption  bands  very  difficult. 

In  the  samples  diffused  at  700°C  for  16  hours  rather  than  the  3  h 
used  above,  the  result  is  the  same  except  the  ^Li  (^Li)  band  at  v  =  350(375)  cm 
has  grown  substantially  with  respect  to  the  other  LVM  bands.  The  other 
bands  retain  their  same  relative  strengths  and  show  little,  if  any,  change. 

In  some  cases  the  changes  were  difficult  to  measure  accurately  because 
relatively  weak  bands  were  close  to  strong  ones  or  nearly  coincided  with 
strong  two-phonon  peaks. 

When  samples  were  diffused  with  /Li  at  800°C  for  one  or  two  hours 
the  350  cm  *  band  increased  dramatically.  A  very  strong  and  broad  band 
at  v  between  380  and  390  cm~^  also  appears  (see  figures  2  (e)  ant*  2  (f)). 

The  rest  of  the  bands  at  V  *  420,  412  and  343  cnT^  show  either  some  small 
growth  or  no  change.  The  ^Li  diffused  samples  also  show  a  large  increase 
of  the  corresponding  band  at  V  *  375  cm“^  and  the  presence  of  the  broad 
band  at  v  between  400  and  410  cm  If  we  take  the  centers  of  the  broad 
bands  at  V  *  381  and  402  cm  \  respectively,  as  a  measure  of  the  band  frequen¬ 
cies,  then  the  small  isotope  frequency  ratio  of  M..055  and  the  broad  band 
width  indicate  that  it  is  due  to  a  new  Li-defect  complex.  It  is  also  indicated 
that  a  significant  part  of  the  kinetic  energy  of  the  vibrational  motion  is 
not  associated  with  the  Li-species. 


Natural  H-dQped  ZnSe  (10Q  ppm)  in  which  Li  was  introduced  during 
crystal  growth  was  also  measured.  It  shows  all  of  the  \i  bands  present 
jtn  the  Li-diffused  samples  except  the  one  at  v  *  420  cm  ^  which  is  too 
weak  to  distinguish  from  the  strong  two  phonon  band  at  v  -  421  cm 
Lithium  was  also  diffused  into  a  CVD  sample  grown  at  Raytheon  Corp.  A 
lower  diffusion  temperature  (600°C)  and  a  longer  time  (24  hours)  was  used 
to  avoid  damaging  the  highly  polycrystalline  sample.  The  result  was  very 
similar  to  that  discussed  above  for  the  grown  material  which  was  diffused 
at  700°C  for  four  hours. 

Lithium  isotope  mixtures  were  used  as  the  diffusion  source  in  an  at¬ 
tempt  to  determine  if  any  of  the  Li  modes  were  coupled  vibrations  involv¬ 
ing  the  motion  of  two  or  more  Li  atoms.  No  new  LVM  bands  were  observed. 
This  indicates  that  all  the  localized  vibration  modes  which  are  observed 
likely  originate  from  defects  involving  only  one  Li  atom/defect.  This 
result  is  completely  different  from  that  for  Li-diffused,  undoped  GaAs, 
where  lithium  forms  at  least  five  different  types  of  defect  complexes. 

All  of  these  complexes  contain  more  than  one  Li  atom  and  the  LVM  are  all 
strongly  coupled  vibrations.  We  therefore  conclude  that: 

(1)  Li  diffused  ZnSe  (undoped)  has  6  bands  related  to  Li  defects 
and  all  bands  show  an  isotopic  shift  indicating  that  the 
vibrational  energy  resides  almost  totally  in  the  Li  atoms. 

(2)  Mixed  Li  isotopic  doping  gives  no  new  bands  indicating  either 
_  one  Li/defect  or  almost  no  coupling  if  there  is  more  than  one 

Li/defect, 

(3)  Li  diffusions  at  different  temperatures  and  for  different 
times  indicate  that  one  of  the  bands  is  not  from  the  same 
defect (s)  as  the  other  five  bands. 


(4)  None  of  the  hands  ohaerved  here  are  the  same  as  those  re¬ 


ported  by  Mitsuishi  et  al,  Moreover  ve  have  never  observed 
the  previous  reported  bands  in  any  of  the  samples  measured  in 
this  study. 

From  these  results  one  can  speculate  concerning  the  nature  of  the  Li  de¬ 
fects  as  follows: 

(a)  Since  the  350(375)  cm  ^  band  appears  to  be  unrelated  to  any 
others  it  must  arise  from  a  Li  defect  with  a  symmetry  which 

is  not  lower  than  tetrahedral.  One  possible  defect  would  be  Li 
which  is  expected  to  be  an  acceptor  with  an  ionization  energy 
of  M).114  eV.  If  this  is  the  case  then  for  the  prolonged 
and  high  temperature  anneals  in  which  the  350(373)  cm  ^  band 
increases  dramatically,  a  similar  increase  in  donor  concentra¬ 
tion  would  be  required.  There  is  however  no  independent  evi¬ 
dence  which  supports  this  possibility.  Moreover,  experiments 

involving  ZnSirAl  +  Li  would  indicate  the  frequency  of  Li  is 

Zn 

considerably  larger  than  350(375)  cm  \ 

(b)  The  remaining  five  bands,  which  show  a  nearly  full  Li  isotope 
shift,  require  either  a  defect  with  two  weakly  coupled  Li  atoms 
per  defect  with  the  Li  atoms  in  inequivalent  sites  or  at  least 
two  differect  types  of  defects  with  one  or  more  Li  atoms  per 


defect  for  each. 


C,  Li  diffused  n^type  ZnSe  (Al.  G a} 

It  is  helieyed  that  Al(Ga)  is  a  substitutional,  shallow  state,  effec¬ 
tive  mass  donor  when  present  in  the  form  of  Al_  .  However,  the  fact  that 

Zn 

the  samples  in  the  as-grown  state  are  highly  compensated,  the  analysis 
of  past  photoluminescence  studies,  and  the  observation  of  six  Al— related 
LVM  bands  in  ZnSe:Al  all  tend  to  indicate  that  Al  participates  in  one  or 
more  defect  complexes.  The  previous  LVM  measurements*'  of  ZnSerAl  and 
ZnSe:Al  counterdoped  with  Cu,  Ag,  and  Au  indicated  that  one  band,  located 
at  v  «  359  cm  ,  had  a  different  behavior  than  the  remaining  bands,  and 
it  was  attributed  to  the  Al^  defect.  The  remaining  bands  are  related 
to  a  complex(s)  involving  Al.  The  annealing  study**  under  controlled 
Se  pressure  discussed  in  part  Ill-a  is  in  agreement  with  this  identifi¬ 
cation  of  the  359  cm  ^  band,  but  this  study  as  well  as  the  previous  ones 
leaves  the  five  remaining _A1- related  bands  in  ZnSerAl  essentially  unex- 
plained. 

There  has  been  a  previous  study  by  Ibuki  et  al?  of  ZnSerAl  +  Li. 
This  study  was  largely  prompted  by  earlier  successes  in  using  Li  as  a 
compensating  species  in  GaAs  which  was  doped  with  a  variety  of  impurities. 
The  Li  compensates  the  crystals  by  forming  complexes,  and  the  Li-related 
LVM  helped  identify  and  characterize  the  defect  complexes.  It  is  known 
from  photoluminescence  studies  of  ZnSe  that  Li^n  is  a  shallow  acceptor  and 
pairs  with  donors  such  as  Al  and  Ga  as  Al(Ga)„  -Li_  .  In  the  Ibuki  et  al. 
study  the  Li  and  Al  were  introduced  during  ZnSe  crystal  growth,  and  the 
spectra  obtained  depended  upon  the  doping  concentrations.  In  their  at¬ 
tempts  to  assign  bands  to  vibrational  modes  of  proposed  defect  centers 
they  were  forced  to  make  some  assumptions  without  adequate  justification. 
However  as  will  be  shown,  most  of  their  assumptions  were  correct.  For 


example,  when  the  Li  concentration  «  Al  concentration,  a  new  band  at 
412(383}  cnT*  for  6Li(.7Li)  was  assigned  to  an  unassociated  Li  acceptor, 
presumeably  Li^,  When  the  Li  concentration  A1  concentration  a  band 
at  exactly  the  same  frequency  was  assigned  to  one  of  the  Li  modes  of  a 
Al„  "Li.  pair.  Moreover,  the  strength  of  the  383  cm  ^  relative  to  the 
other  bands  of  the  pair  was  different  from  that  for  the  412  cm  *  in  the 
high  Li  concentration  case  making  the  isotopic  identification  question¬ 
able.  The  same  type  of  situation  existed  for  the  359  cm  ^  band  which  at 
high  Li  concentration  became  assigned  to  an  A1  vibration  of  Al^-Li^  Palrs 
rather  than  the  triply  degenerate  vibration  of  isolated  AlZn  donors. 

The  above  results  indicated  that  a  further  study  of  ZnSe:Al  +  Li 
is  desireable.  We  have  made  a  number  of  measurements  for  this  system  and 
for  ZnSe:Ga  +  Li  as  well.  _  Illustrative  spectra  for  each  of  these  cases 

are  shown  in  Fig.  3  at  the  end  of  this  section.  The  samples  were  Li-diffused  ZnSe 
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doped  with  300  ppm  of  A1  or  an  unknown  concentration  of  Ga.  The  experi¬ 
mental  results  and  qualitative  conclusions  can  be  summarized  as  follows: 

(1)  ZnSe:Al  samples  which  were  annealed  under  different  Se  pressures 

(see  part  III-A  of  this  report)  prior  to  Li  diffusion  showed 

bands  after  Li  diffusion  which  were  reduced  in  the  same 

manner  as  the  5  Al  bands  (not  the  Al_  at  359  cm  )  in  the 

Zn 

pre-diffused  sample.  This  includes  the  post-diffusion  359  cm 
band  thus  indicating  that  this  band  may  indeed  be  due  to 

The  figure  from  the  supplier  was  3000  ppm  however  the  intensity  of  the 
observed  LVM  spectra  would  seem  to  indicate  a  much  lower  concentration. 


ABSORPTION  COEFFICIENT 


different  defects  in  ZnSe;Al  and  ZnSe;Al  +  Li. 

(21  The  352  and  35Q  hands  which  appear  in  both  ZnSe;Al’-^Li(^Li) 
samples  do  not  appear  in  either  ZnSe;Ga-\i(^Li)  and  thus  these 
hands  are  attributed  to  A1  modes  of  the  AlZn-Li2n  P3***  Thi3 
conclusion  is  in  agreement  with  Ibuki  et  al. 

(3)  Six  bands  are  observed  in  the  Al  +  ^Li  samples  and  only  five 
in  Al  +  \i  samples.  With  the  exception  of  the  359  cm  ^  dis¬ 
cussed  above  all  of  the  bands  differ  in  frequency  from  those 
in  the  Al-doped  samples.  Here,  as  in  the  Ibuki  work,  the 
383  cm  band  of  the  Al  +  ^Li  spectrum  was  much  larger  than 
expected  by  its  isotopic  counterpart  at  412  cm  ^  in  Al  +  ^Li 
samples.  This  suggests  the  reason  that  only  five  bands  are 
observed  in  the  Al  +  ^Li  spectrum  is  due  to  the  383  cm  ^ 
absorption  band  is  the  result  of  the  superposition  of  two  bands, 
one  of  which  is  the  412  cm  1  band  in  the  ^Li  diffused  samples. 

(4)  In  comparing  the  Al  +  ®Li  and  Al  +  \i  spectra,  two  bands  show 

no  shift  (350  cm  ^  and  359  cm  and  two  show  almost  a  full  Li 
isotopic  shift  (420  cm  1  -*■  392  cm  ^  and  412  cm  ^  383  cm  ^).  The  remain¬ 

ing  two  bands  show  substantially  less  than  the  full  isotopic 

shift  (389  cm  \  373  cm  ^  383  cm  ^ ,  356  cm  indicating 

compled  modes  involving  both  Al  and  Li  vibrational  motion.  This  result 

suggests  the  Al^-Li^  type  of  defect  proposed  in  photoluminescence 

studies  and  a  model  based  on  this  defect  is  applied  to  both 

the  Al_  -Li_  and  Ga  -Li  observations  in  part  IV  of  this  report. 


D.  Hardness  of  ZnSe;  Annealing  Effect 

As  mentioned  previously  the  samples  purchased  from  Eagle-Picher 

Corp.  are  polycrystalline  with  large  size  grains,  The  grain  size,  which 
*3  -3 

is  approximately  10  volume  £  10  mm  ,  depends  somewhat  upon  the  dopant 
concentration.  Indentations  were  made  on  the  central  area  of  each  grain 
in  order  to  minimize  grain  boundary  effects.  Measurements  were  made  with 
undoped  and  doped  samples  to  determine  if  there  is  an  impurity  effect. 

1*  Undoped  ZnSe 

The  procedure  employed  is  first  to  make  at  least  five  groups  of 
approximately  five  indentations  each  aligned  parallel  in  an  arbitrary 
direction,  then  rotate  the  sample  90°  and  make  another  similar  five 
groups  of  indentations.  The  average  of  more  than  twenty  readings  in  one 
direction  was  computed.  The  result  based  on  five  different  large  size 
grains  from  two  ingots  indicates  that  within  the  experimental  accuracy 
the  material  is  isotropic.  The  K.H.N.  is  77.2  and  the  scatter  of  the 
readings  is  less  than  ±  4%.  In  regions  close  to  the  grain  boundaries, 
the  hardness  is  greater  by  about  15%  20%. 

2.  Al~doped  ZnSe;  100  ppm  and  300  ppm 

For  heavily  Al-doped  ZnSe,  chemical  polishing  reveals  some  "parallel 
strips"  on  the  surface.  The  origin  of  these  strips  is  not  known,  and  we 
have  not  had  the  opportunity  to  investigate  them.  Indentations  were 
made  along  the  directions  perpendicular  to  (X)  and  parallel  to  (If )  the 
stripes.  In  all  but  one  case,  the  K.H.N.  (X)  is  smaller  than  the  K.H.N. 
(II).  For  the  more  heavily  doped  ZnSe:Al  (300ppm),  the  difference  in  the 
K.H.N.  between  the  two  directions  is  about  10%  while  the  scatter  of  the 
measurements  for  a  given  condition  is  generally  with  ±  4%.  The  K.H.N. 


values  are  shown  In  Table  1,  The  results  are  compatible  with  previous 
measurements  we  made  which  used  heavier  loads,  e.g.  load  >  200  gm.  The  results 
indicate  that  ZnSe  is  hardened  by  the  addition  of  A1  and  the  K.H.N.  in¬ 
crease  can  be  by  as  much  as  4Q%  for  ZnSe:Al  with  an  A1  concentration 
-  300  ppm. 

3.  Annealing  effect  ZnSe:Al  (300  ppm) 

Hardness  measurements  were  made  on  the  heat-treated  samples  to  deter¬ 
mine  the  effect  of  annealing  on  the  microhardness.  Usually  the  samples 
have  some  slight  surface  deterioration  after  annealing,  and  need  repolish¬ 
ing.  There  is  no  measureable  difference  in  the  K.H.N.  value  before  and 
after  repolishing  the  samples  surface.  Results  from  three  different  in¬ 
gots  (with  A1  300  ppm)  are  shown  in  Table  2.  The  microhardness  in 
both  directions  decreases  after  annealing  at  900°C  in  a  Se  partial 
pressure  of  2  atms.  After  the  first  40  hours  of  annealing,  the  K.H.N. 
values  decrease  about  10%,  and  the  amount  of  decrease  in  the  next  40  hours 
of  annealing  is  less  than  5%.  There  is  a  similarity  between  the  observed 
decrease  in  the  microhardness  and  the  decrease  in  the  band  strength  of 
five  unknown  Al-related  bands.  It  indicates  that  there  may  be  a  simple 
relationship  between  the  defects  associated  with  five  unknown  Al-related 
bands  and  the  defects  associated  with  the  change  in  the  microhardness. 

There  is  however,  insufficient  data  to  establish  a  quantitative  correla¬ 
tion.  As  previous  mentioned  the  LVM  spectra  indicates  the  presence  of  a 
defect  complex(s)  involving  Al  but  the  specific  complex(s)  could  not  be 
Identified  with  certainty.  If  the  complex  were  known  then  we  are  also 
unaware  of  any  theory  relating  the  microhardness  of  semiconductors  to  the 
concentrations  of  specific  types  of  defect  species.  Therefore  only  the 


"Table. 


following  qualitative  conclusions  can  he  drawn; 

1  The  addition  of  A1  hardens  ZnSe 

2  Annealing  Al-^doped  ZnSe  can  reduce  the  microhardness 

3  The  reduction  in  microhardneaa  in  the  annealing  may  be  related 
to  the  changes  in  an  Al-related  defect  complex(s)  concentration. 


IV.  Model 


In  this  section  a  simple  yalence^farce  type  of  model  calculation 

will  be  given  for  the  Al_  -Li_  and  Ga_  ~H_  defects*  First  let  us 

z,n  zn  Zn  Zn 

briefly  review  the  experimental  features  to  be  fitted  by  the  model.  The 

6  7t  . 

A1  +  Li  doped  ZnSe  has  6  absorption  bands  and  the  A1  +  L1  has  only  5 

bands.  Comparing  the  two  spectra  it  appears  that  two  modes  do  not  change 

frequency  at  all,  the  two  highest  frequency  modes  for  ^Li  show  almost  a 

"full  Li  isotopic  shift"  when  ^Li  is  replaced  by  \i,  and  the  remaining 

two  modes  show  only  a  partial  shift  with  one  mode  superimposing  on  one 

of  the  modes.  We  therefore  anticipate  that  two  modes  will  be  Al-like, 

two  Li-like,  and  two  will  be  coupled  modes.  If  the  A1  is  replaced  by 

Ga  in  the  model  the  larger  Ga  mass  can  be  expected  to  greatly  reduce  the 

A1  mode  frequencies  so  that  they  merge  with  the  band  mode  continuum  and 

are  no  longer  LVM.  Thus  one  expects  only  three  Li-like  modes  with 

frequencies  which  show  almost  a  full  isotopic  shift. 

The  model  used  to  estimate  LVM  frequencies  for  ZnSe:Al(or  Ga)  +  Li 

employs  the  following  assumptions: 

a)  The  A1  is  a  substitutional  donor  on  a  Zn  site,  Al_  ,  and  Li  acts  a 

Zn 

compensating  acceptor  LI  a  2—  nearest  neighbor  site. 

Zn 

b)  Since  the  masses  of  Li  and  A1  are  much  less  than  those  of  Zn  and  Se, 
we.  treat  both  impurity  atoms  as  vibrating  against  a  rigid  lattice. 

c)  The  forces  acting  on  the  Al(Li)  are  between  it  and  its  four  tetra- 
hedrally  arranged  Se  nearest  neighbors  plus  the  interaction  between 
the  A1  and  Li.  For  this  simple  model  we  neglect  forces  from  more 
distant  neighbors. 

d)  The  interaction  between  the  A1  and  Li  is  the  sum  of  the  attractive 
Coulomb  force  and  a  short  range  force.  The  force  constant  is  treated 
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as  an  experimentally  determined  factor,  and 
e)  Lattice  distortion  is  introduced  by  changing  the  bond  angles  between 
the  Aiai)  and  its  four  surrounding  Se  atoms.  The  distortion  reduces 
the  symmetry  of  the  potential  at  the  Al(Li)  site  and  produces  the 
frequency  splitting.  The  small  but  non-zero  isotope  shift  for  two 
of  the  modes  (see  part  III-C  of  this  report)  indicates  that  a  strong 
coupling  model  with  the  changes  in  the  equilibrium  positions  of  the 
A1  and  Li  should  also  be  introduced. 

Fig.  4  defines  the  coordinate  system.  The  angles  and  atomic  coordi¬ 
nates  are  also  given  in  the  figure.  The  black  atom  #1  at  (0,0,0)  reprte- 
sents  the  A1  and  #2  is  the  Li.  The  "a"  is  the  lattice  constant  of  ZnSe. 
The  white  atoms  #3  -*■  #9  represent  the  Se  atoms  which  are  the  nearest 
neighbors  to  the  impurities. 

The  equations  qf  motion  are 


IM-  r  -  p-r 


where  i  labels  the  atoms,  "1"  for  the  A1  and  "2"  for  the  Li,  and  a  labels 
the  Cartesian  coordinates.  The  m^  is  the  mass  of  the  i  th  atom,  ria  is 
the  a:  component  of  the  displacement  from  equilibrium  of  the  i  th  atom. 
F^g  represents  the  coupling  constant  between  the  a  component  of  the  dis¬ 
placement  of  the  i  th  atom  and  the  8  component  of  the  j  th  atom,  and  is 
defined  as 


_  __L 

A 


where  <P(i,j)  is  the  interaction  potential  between  the  particle  i  and  j. 
jP  is  a  symmetrical  6x6  matrix  which  can  be  separated  into  four  submatricies. 
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where  (F^g) ,  (F*g),  (F^g)  and  (F^g),  are  four  3x3  matricies. 
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The  lower  left  hand  submatrix  F^g  is  the  transverse  of  the  upper  right 


.12 


hand  quarter  F~~.  Both  of  them  represent  the  coupling  force  matricies 


between  particles  1  and  2.  We  assume  that  the  interaction  between  Al  and 

Li  is  the  sum  of  the  attractive  Coulomb  potential  <pc  and  a  short  range 

,12 
ae 


s  12 

potential  <Jr.  Then  the  elements  of  the  upper  right  hand  quarter  FnR  can 


be  written  as 


F  A  *  -  - - -  (  <t>e  4  d>s) 

9r(<  9  rtfl  *  J 


r*r„  j 


C  g 

Since  both  $  and  4>  are  functions  of  the  distance  between  particles 


1  and  2,  on^y»  i»e*  $  -  $  (r,,)  and  4>S  *  $S(r10).  Then  we  can 
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k12' 


calculate  the  derivatives 
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and 

3*<t>c  j  Q 

3X|  *  »*«  3Z2  "  d  X,  5Yi  ~  3Xi  3Z2 

Where  and  *  the  "effective  charge"  of  the  A1  and  Li,  €  *  the 
relative  dielectric  constant,  rQ  *  the  equilibrium  distance  hetween  the 

SU  9t  3 

A1  and  Li,  <f  and  $  ■  the  second  and  first  derivatives  of  9  with 

respect  to  the  interatomic  distance  We  can  see  that  only  the  dia¬ 

gonal  terms  are  net  equal  to  zero.  So 


T^ie  Fxx»  Fyy*  an<*  Fzz  be  determined  by  experiment. 


The  upper  left  quarter  JF^g  which  represents  the  restoring  force  ma¬ 
trix  can  be  split  into  the  sum  of  Al^e  forces  and  Al^i  force, 

Fc*g  *  (  VA1-Se  +  ^ot^Al-Li1 

11  12 

where  R  "(F^g) ,  We  assume  that  Al-Se  forces  are  hamonic 

forces  with  force  constant  K^,  then 


(F^W=  4-^0), 

J  =  3,415,6  ,K  JP 


where  r  *  the  interatomic  distance  between  particles  1  and  j  and  the 
summation  is  over  four  surrounding  Se  atoms.  For  example,  the  derivative 
can  be  calculated,  and  then  we  get 


(F  ii  )AI-S%  -  "K>  to*ae  -  k,  cosae,  , 

(Fm)am.  =  °‘, 


(  F  i  £  )ai-5€  =  -  K,  cos  e  Sin  e  -  ic,  Sin  0|  Cos  6, 


The  equation  of  motion  of  particle  1  in  the  x  direction  is  then 


Wi  Xi  =  21  FxJ  F j£ 

jalyi 


-  £  Ki  ( c«ae  +  t«je, )  +  &"  + 

-  £  K|  +  Jc ,  cos 0j sinfijJ  2,  ”h  (  "F  4>s 
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Lattice  distortion  is  introduced  hy  lettine  the  ai  c 

V  J-etUng  the  Al~Se  and  Li-Se 
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OTd  «■-«•<*»  bond.,  The  elMMt8  „f  cll.  , 

e"e”t8  of  che  f"«  constant  tan  be  written 

as  follows : 
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where  *  the  force  constant  of  the  Al-Se  bond 

«  the  change  in  the  force  constant  of  the  Al-$e  bond 
^  ■*  the  force  constant  of  the  Lt-Se  bond 
AK^  ■  the  change  in  the  force  constant  of  the  Li-Se  bond 

Because  the  Cs  symmetry,  the  vibrations  perpendicular  to  the  mirror 

plane  are  uncoupled  from  the  vibrations  in  the  plane.  The  eigenfrequencies 

cies  and  eigenvectors  are  obtained  by  diagonalizing  a  2  x  2  matrix  and 
2 

a  4  x  4  matrix  B  ,  respectively,  where 
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If  we  assume  that  when  there  is  no  Al-Li  coupling  and  no  lattice 

distortion,  then  the  two  triply  degenerate  eigenfrequencies  are  the  local 

mode  frequencies  of  Al^n  and  Li^,  respectively,  which  then  yield  and 

as  experimentally  determinded  parameters.  With  no  Al-Li  coupling 

the  LVM  frequencies  are  taken  to  be  359  cm  *  for  Al_  (T,  symmetry)  and 

zn  a 

411  cm  *_for  ^Li„  (T .symmetry) ,  These  choices  for  the  uncoupled  fre- 
Zn  a 

quencies  were  discussed  in  earlier  sections  of  this  report. 

When  the  Al-Li  coupling  is  introduced  the  calculation  proceeded  as 
follows : 


From  equations  4-3  and  4-12  we  note  that  the  simple  lattice  distortion 


did  not  Alter  ttxe  force  constants  of  AJL^  and  Liy,  where  denotes  the 
ft  component  of  A  Yihration,  If  we  further  assume  <5  «  0,  then  the  A1 


and  Liy  are  decoupled  from  each  other  and  would  have  the  same  frequencies 
as  v(Alz^)  and  vCLtZn) ,  respectively.  The  359  cm'"'*'  band  is  thus  assigned 
to  the  local  model  of  Al^  and  the  411  (or  383)  cnT*"  band  to  the  local  mode 


2 

The  A  x  A  matrix  B  is  diagonalized  to  get  the  eigenfrequencies  and 
eigenvectors  of  the  coupled  modes  in  the  mirror  plane.  The  local  mode 
frequencies  were  computed  as  a  function  the  x  component  of  the  Al-Li 
coupling  constant  y  in  fig,  5.,  The  paramenters  9,4>,  AK^  and  Al^  are  cho¬ 
sen  so  that  when  y  *  0,  the  frequencies  follow  the  order,  v(Li  )  >  v(Li  ) 

v(Al  )  >  v(Al  ),  Because  with  6^0  and  6  ■  0  the  Al-Li  coupling  is 
X  z 

only  in  the  x  direction,  modes  1  and  A  remain  essentially  Li^-like  and 

A1  -like,  respectively.  Their  eigenfrequencies  decrease  slightly  as 
z 

Y  increases  while  chose  for  modes  2  and  3  decrease  rapidly.  When  the 
coupling  is  weak,  v  (mode  2)  decreases  faster  than  v  (mode  3)  because  the 
former  consists  primarily  of  Li^  vibration  and  the  latter  primarily  of 
Alx  vibration.  But  as  the  coupling  increases  and  v  (mode  2)  approaches 
v  (mode  3)  the  two  modes  become  strongly  mixed  and  there  is  no  cross¬ 
over  of  the  frequencies.  -When  the  coupling  is  weak  the  character  of  the 

mode  2  and  3  is  that  mode  2  is  Li  -like  and  the  mode  3  is  A1  -like.  For 

x  x 

sufficiently  strong  coupling  (Y  large)  they  reverse  character.  Figure  6 

shows  the  kinetic  energy  composition  of  the  coupled  modes  2  and  3, 

KE(Li) 

i.e.  KE(to~tai)  *  as  3  function  of  the  Al-Li  coupling  constant  Y. 

If  we  replace  the  mass  of  A1  by  Ga  while  keeping  the  other  parameters 
unchanged,  then  we  will  get  the  coupled  modes  for  the  Ga-Li  pair. 

Because  of  the  larger  mass  of  Ga,  the  lower  two  modes  which  consist 
primarily  of  Ga  vibration  merge  in  the  fundamental  absorption  band  and 
do  not  appear  in  Fig.  5.  The  dashed  line  in  the  same  figure  shows  the 
higher  two  modes  1’  and  2'  which  consist  primarily  of  Li  vibration. 
Because  the  large  energy  difference  between  the  Lix  and  Gax  vibration, 
mode  2'  remains  Lix~like  thr  jghout  the  y  range  in  this  figure. 


WAVE  NUMBER  (cm-l) 


We  can  estimate  the  Al-Li  interaction  from  Fig.  5.  The  parameters  for 

the  best  fitting  of  experimental  data  are  given  in  Table  3.  The  amount  of  the 

lattice  distortion  in  the  Li  site  is  less  than  the  A1  site.  This  is  compatible 

o 

with  the  fact  that  the  tetrahedral  covalent  radius  of  A1  is  1.26  A  and  the 

o 

ionic  radius  of  Li  is  2.26  A.  Thus  the  attractive  Coulomb  interaction  between 
die  A1  and  the  Li  could  procude  a  new  A1  equilibrium  position  which  is  further 
from  the  original  T^  site  than  that  for  the  Li. 

In  the  Al-^Li  system,  we  only  replace  the  mass  of  ^Li  by  ^Li  while  keeping 
all  the  parameters  unchanged.  The  frequency  of  mode  2"(^Lix-like)  is  smaller 
than  that  of  the  mode  3"  (Al^-like)  even  when  the  Al-Li  interaction  y  is  zero 

as  shown  in  Fig.  7.  As  we  increase  the  interaction,  the  frequency  of  mode  2" 
decrease  with  a  faster  rate  than  the  mode  3".  Again  the  model  frequencies 
do  not  cToss-over,  however  in  this  case  when  there  is  strong  Al-Li  coupling, 

mode  3"  remains  A1  -like  and  the  mode  2"  is  Li  -like.  Figure  8  shows  the 

X  x 

KE(Li) 

kinetic  energy  composition  of  the  coupled  mode  2"  and  3",  i.e.  » 

as  a  function  of  the  Al-Li  coupling  constant  y.  If  Y  becomes  sufficiently 
large,  there  will  be  a  cross-over  region  for  the  modes  2"  and  4"  (Alz~like) . 

The  mode  2"  will  change  from  Lix~like  to  Alz-like  and  mode  4"  will  go  from  Al^- 
like  to  Lix~like. 

This  simple  model  explains  that  the  389  and  373  cm  ^  bands  in  ZeSe:Al  +  Li 
are  the  strong  coupling  modes  while  their  corresponding  383  and  356  cm  ^  bands 
in  ZeSe:Al  +  \i  remain  as  A1  -like  and  Li  -like  modes,  respectively. 

X  x 


T«Wle 


WAVE  NUMBER  (CM-I) 


Y  X  I04  (  cm-2  x  No.  of  H  atoms  ) 
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